Addendum: Behavior of a bipartite system in a cavity 
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This note is an Addendum to our previous article [Phys. Rev. A 81, 053820 (2010)]. We show that 
under the assumption of a Bose-Einstein distribution for the thermal reservoir, zero-temperature 
properties of the entangled states considered there are not changed by heating, for temperatures up 
to the order of room temperatures. In this case, the system is dissipative in free space and presents 
stability for a small cavity, both for T = and for finite temperature. 
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Introduction: In this note we consider how heating 
could affect the properties of the bipartite system studied 
in [l| . This bipartite system consists of two dressed atoms 
inside a reflecting cavity which do not interact with each 
other, only with an environmental field. Our approach to 
this problem makes use of the notion of dressed thermal 
states 0, in the context of a model already employed 
in the literature, of atoms, or more generally material 
particles, in the harmonic approximation, coupled to an 
environment modeled by an infinite set of pointlike har- 
monic oscillators (the field modes). The dressed ther- 
mal state approach is an extension of the dressed (zero- 
temperature) formalism which has been introduced ear- 
lier. The reader is referred to [l], [|J for details of our 
formalism. 

We here consider entanglement as a pure quantum ef- 
fect, a characteristic of quantum mechanics, which is also 
nonlocal, in the sense that distant and nonintcracting 
systems may be entangled. This is due to the physical 
meaning attributed to superposed states, a concept with 
no correspondence in classical physics, and not to the 
interaction between the (in our case, dressed) atoms. In- 
deed, such properties of entanglement of noninteracting 
systems have been used in the realm of teleportation and 
quantum information theory and to conceive quantum 
communication devices. We think that the investigation 
of in which measure the properties of such systems are 
affected by heating is important, particularly in the case 
of room temperatures. The possibility of constructing 



such devices at temperatures of everyday life would be 
an interesting matter. 

In this note we remain on theoretical grounds, and 
perform a study of thermal effects on our bipartite sys- 
tem, assuming that the dressing fields of the atoms obey 
the Bose-Einstein distribution and that the first-level ex- 
cited dressed atomic states evolve in the same way as 
in [l|. Under these assumptions, it results that the re- 
duced density matrix does not depend on temperature. 
On the other hand, employing methods and relying on 
results obtained in 0, y|, one finds that the individual 
dressed atoms in the bipartite system are very little af- 
fected by heating for temperatures around 300 K. The 
overall conclusion is that, methods and results from [l[ 
are valid if we consider the system at room temperatures. 

Our bipartite system is composed of two subsystems, 
A and B; the subsystems consist respectively of dressed 
atoms A and B and the whole system is contained in a 
perfectly reflecting sphere of radius R in thermal equilib- 
rium with an environment (field) at a temperature /3 _1 . 
We consider each atom carrying its own dressing field 
cloud, independently of each other. This means that we 
are considering the behavior of a noninteracting bipar- 
tite system. Also, we take an specific type of reservoir, 
in the thermal Bose-Einstein state, where the excitation 
distribution is diagonal in the number basis. 

The density matrix: We start, at time t = 0, from a 
family of superposed states of the bipartite system given 
by 
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where < £ < 1 
In Eq. 0, 
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\0a1b) stand respectively for the t = states in which 
the dressed atom A (B) is at the first level, the dressed 
atom B (A) in the ground state. They are 
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(o)) = |rf(o))®|r B (o)) 
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and 
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(O)) = |r 4 (o))®|rf(o)). 
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In order to take into account temperature effects, we 
consider the density matrix, 

PAis(t, 0) = \%ab) (*ab| <8> pa{P) ® Pb(P), 

in terms of dressed objects. We take the density matrix 
for the thermal bath obeying the Bose-Einstein distribu- 
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tion, 

pAP) = T\^re- h ^ a l a * 
\ Zk 

1 oo 

= 11^ E e- H ^ n "\n k )(n k \, (4) 

k n k —0 

where n k = 0, 1, 2, . . . correspond to the occupation num- 
I 



where we have changed notation, |r^(i)) = and 
similarly for the other states. 

Putting together Eqs. © and ([7]), the reduced den- 



bers of the dressed field modes fc, and Z k is obtained 
imposing the condition of unit trace for pyi(/3); this gives 



and so, 
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sity matrix is obtained by taking the trace over the field 
modes, 



g — ft/3(wini+U2«2+---) 

= E ZiZz |ni)(m|® |n 2 )(n 2 |(8)---. 

ni,ri2,... 

For the part containing only the atoms, we have 

\*ab)(Vab\ = €[|U(t))(lA(t)|®|0 B )<0 fl |] 

+(1 -0 [10a ) (0^1 <8 |l B (t)> <1 B (*)|] 

+ V^W)e l * [|0^> ® \l B (t)) (0 B |] 

+ VW^T)e-^ [|1a(*)> (0a| ® |0 B ) (l B (t)|] , 
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To calculate the matrix elements above we use the time A similar equation holds for |1#(£); m, 112, ■ ■ ■ )• Inserting 

Eq. © and its equivalent for the atom B into Eq. [5J , we 
obtain, after some rather long manipulations, 



evolution of the states involved [1] 

\l A (t);n 1 ,n 2 ,...) = ^/^(i) |1„(0); n u n 2 , ■■ ■) , (9) 



where the coefficients fAv(t) are 

fA„{t)=Y,%tle- iU " t ; £l/^(*)| 2 = i 



(10) 




\fAA(t)f $p A 0$r A 



3p B odr B 



i=l 

+ V 7 C (1 - l)e^/lA(*)/BB(*)^O^rAl^fll 5 r B + V 7 ^ 1 - O e_ ^/AA(i)/BB(i)^Al^ A 0<5pBO^ B l} 



But 



^7, 



= 1, 



v», 



(11) 



and, therefore, the temperature-dependent terms disap- 
pear, leaving the reduced density matrix with exactly 
the same form as in For identical atoms we have 



f AA {t) = f BB (t) 


= foo(t) and the nonvanishing matrix 


elements of p are 




pM 


= 1-l/ooWI 2 , 


Pol(t) 


= (l-0l/oo(*)| 2 , 


Ploit) 


- el/ooWl 2 , (12) 




= Ve(l-e)e^|/oo(*)| 2 , 




= \/£(l-Oe- i0 |/oo(«)| 2 . 



Notice that this property is not specific to the Bose- 
Einstein distribution we have employed, but it is com- 
mon to any normalized distribution which is diagonal 
in the number basis. This may not happen for other 
types of reservoir, like coherent squeezed thermal states, 
where the factorization of the temperature-dependent 
terms may not occur. 

Conclusions: At finite temperature, all the results that 
depend only on the reduced density matrix remain the 
same as in the zero-temperature case This conclu- 
sion is valid for our choice of the Bose-Einstein distri- 
bution to the heated environment, using Eq. © to the 
time evolution of the first-level excited atomic states. 
The choice of the Bose-Einstein distribution to the field 
modes can be justified: in the case of an arbitrarily large 
cavity, the dressed field modes coincide with the bare 
ones [2] and in the limit of vanishing coupling with the 
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atom, these modes follow the Bose-Einstein distribution 
exactly. Strictly speaking, this is not the case for the 
coupled atom-field system in a finite cavity. Neverthe- 
less, in many situations this approximation is acceptable 
in the weak-coupling regime [4{. The use of Eq. (|9]) to 
the time evolution of the first-level excited atomic states, 
as we will see below, is a good approximation for room 
temperatures. 

In this case, heating does affect each atom individually. 
In Q and in 0], some of us have studied the time evo- 
lution of the temperature-dependent occupation number 
of a single dressed atom, n' Q (t,/3), for respectively a very 
large cavity (free space) and confined in a small cavity, 
with the same heating procedure as we have used here. 
In the case of a large cavity, we have shown that, starting 
from an initial value n' (t = 0,(3) = 1, n' (t,f3) evolves, 
for weak coupling, steadily from 1 to an equilibrium value 
of 

<(oo,/3) = l/(e r ^-l). 

For T = 300 K, for a frequency ui = 4.0 x 10 14 /s, we have 
n' (oo,T — 300) w 0.09, a value slightly higher than the 
corresponding value for T = 0, n' (oo,T = 0) = 0. For 
the same initial condition, n' (t = 0, fi) — 1, and the same 
value of the emission frequency and temperature as in the 
large cavity, it is found in Q for a small cavity (radius 
R ~ 10 -6 m) that the time evolution of the occupation 
number has an oscillating behavior. Raising the tem- 
perature increases the average value of the occupation 
number, but this increase is significative only for high 
laboratory temperatures. For instance, for T = 10 5 K 
(~ 8.4 eV) (a value lower than the ionization tempera- 
ture of 13.6 eV for the hydrogen atom) for a frequency 
lu = 4.0 x 10 14 /s, as above, one obtains that the aver- 
age occupation number is about 3.5 times higher than 
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the zero-temperature value, n' (t, T = 0) ~ 1. On the 
other hand, this effect is negligible for room tempera- 
tures; in this case, taking for instance a temperature 
T > 300 K, it is found that the occupation number is 
given by n' (t, T > 300) > 1, which will remain very close 
to the zero-temperature value. Therefore, if we consider 
the system at room temperatures, individual atoms are 
very little affected by heating and the use of Eq.© is 
justified. 

Our system is composed of two noninteracting atoms, 
each one carrying its own dressing field. The superpo- 
sition principle dictates the existence of correlations be- 
tween the parts of the system (the two atoms) even if 
they do not interact. By introducing temperature in 
the environment by means of the Bose-Einstein distri- 
bution, the result is that the entanglement properties of 



the considered bipartite system are not affected by heat- 
ing up to room temperatures. In particular, considering 
two identical atoms, the concurrence, the entanglement 
of formation and the negativity depend only on the ma- 
trix elements in Eq. tjl 2|) and thus, cannot depend on 
the temperature. However, the temperature of the en- 
vironment cannot be so high that the atoms dissociate, 
or such that the adopted atomic model becomes invalid. 
More precisely, from a physical point of view, our model 
applies for room temperatures. As we have seen above, 
for temperatures of the order of ~ 300 K or not much 
higher, the harmonic approximation we have used in [l[ 
remains valid and the time evolution of the thermal oc- 
cupation number of the atoms is very close to the zero- 
temperature case. This establishes the range of validity 
of the present study. 
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